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We review the current theoretical situation with regard to the anomalously large cross section for 
exclusive J/ip + rjc production in e^e~ annihilation, as measured by BELLE Collaboration. 
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I. INTRODUCTION 

The J/ip, a J^'~' = 1~~ charmonium state, has, for 
many years, been a nice probe of hadron physics. Exper- 
imentally, it is easily detectable through its decay into a 
lepton pair. Theoretically, it has been a testing ground 
of perturbative quantum chromodynamics (QCD) be- 
cause it involves both long and short-distance dynamics. 
Because of the separation between the long and short- 
distance scales involved in quarkonium dynamics, one 
can make use of nonrelativistic QCD (NRQCD) [1] to 
describe the production and decay of heavy quarkonium 
in a factorized form that is analogous with those in stan- 
dard hard-scattering QCD factorization theorems [2-4]. 
The factorized formula is a linear combination of long- 
distance NRQCD matrix elements, which scale according 
to a known power of the velocity v of the quark inside 
a quarkonium in the meson rest frame. The matrix ele- 
ments are determined by experimental measurements or, 
in some cases, through lattice calculations. Once the uni- 
versal long-distance matrix elements are determined, one 
only needs to calculate the corresponding perturbative 
short-distance factors in order to predict the value of a 
physical observable. 

As the first effective field theory for treating heavy 
quarkonium physics, NRQCD succeeded in providing 
infrared-finite predictions of the P— wave quarkonium de- 
cay rate [5] and presented a solution to the excessively 
large production rates for J/ip and ip' at large transverse 
momenta at the Tevatron [6] . These results were enough 
to rule out the color-singlet model [7], which had sur- 
vived for decades. On the other hand, there are open 
questions still to be resolved in NRQCD phenomenology. 
One is the polarization of the J/ijj at large transverse 
momentum at the Tevatron: the CDF collaboration has 
not observed the large transverse polarization [8] at large 
transverse momentum that was expected from theoreti- 
cal predictions [9-11]. However, the experimental data 
are not yet decisive. Experimental analyses with high- 
statistics, Run-II data are under way. On the theoretical 
side, large relativistic corrections have been found [12]. 
Furthermore, predictions have, so far, ignored spin-flip 
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interactions. Lattice calculations to estimate the impor- 
tance of these interactions are being carried out [13]. 

Recently, quarkonium physics has been faced with two 
particularly difficult problems, both of which arise from 
BELLE measurements of J/tp production in e'^e~ anni- 
hilation at y/s = 10.6 GeV [14]. One production rate for 
inclusive J/ip + cc relative to J/ip + X [14], which is much 
larger than predictions [15] . The other is the cross section 
[14] for the exclusive process e+e~ ^ J/ip + rjc, which 
is larger than predictions by at least an order of magni- 
tude [16, 17]. Both of them are the largest discrepancies 
that currently exist in the Standard Model. 

In this proceeding, we review current theoretical status 
with regard to the exclusive problem. We first consider 
the original discrepancy between the data [14] and the 
theory [16] for exclusive J/Tp + rjc production. Then we 
summarize two recently proposed scenarios to reduce the 
discrepancy between the data and the theory. One sce- 
nario is that the BELLE signal may contain double- J/-)/) 
events[18, 19]. The other is a conjecture that the signal 
may include exclusive J/ip production associated with a 
glueball [20]. 



II. EXCLUSIVE J/V + vc PRODUCTION 

According to the NRQCD factorization formalism, ex- 
clusive processes of quarkonium production and decay, in 
which there is no other hadrons, should be described by 
the color-singlet model, up to corrections that are higher 
order in w. As in the simplest examples, such as elec- 
tromagnetic annihilation decays and exclusive electro- 
magnetic production processes, e^e~ annihilation into 
exactly two charmonia should be accurately described 
by the color-singlet model, too. Because of the mono- 
energetic nature of a two-body final state, the absence of 
additional hadrons in the final state can be guaranteed 
experimentally. For many charmonia H, the NRQCD 
matrix element can be determined from the electromag- 
netic annihilation decay rate of either H or of another 
state that is related to H by the heavy-quark spin sym- 
metry. Cross sections for exclusive two-charmonium pro- 
duction in e+e^ annihilation can, therefore, be predicted, 
up to corrections that are suppressed by powers of w^, 
without any unknown phenomenological factors. In con- 
trast inclusive quarkonium cross sections in hadroproduc- 
tion involve several matrix elements whose values have 
large uncertainties. 
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FIG. 2: QED diagrams that contribute to the color-singlet 
process 7* —^ cci(^S'i) + cci. 

Unfortunately, cross sections for exclusive two- 
charnionium production were too small to be measured 
until the newly built high-luminosity B factories collected 
large-statistics data. A naive estimate of the cross sec- 
tion for J/'ip + r]c in units of the cross section for 
is 

Rm+vc] ~ i^Y- (1) 

V .^beam / 

The two powers of are the minimum needed to pro- 
duce a cc + cc final state. There is a factor of (nicv)^ 
associated with the wave function at the origin for each 
charmonium. These factors in the numerator are com- 
pensated by factors of the beam energy Ebeam in the 
denominator, which yield a dimensionless ratio. As 
an example, consider e~^e~ annihilation with center-of- 
mass energy 2£Jbeam = 10.6 GeV. If we set v"^ « 0.3, 
as « 0.2, and nic « 1-4 GeV, we obtain the naive esti- 
mate R[J/t/j + ?7c] ~ 4 X 10^^. This should be compared 
with the total ratio it![hadrons] ~ 3.6 for all hadronic 
final states [21]. 

Because the two charmonia in the final state have op- 
posite charge-conjugation quantum numbers, it is as- 



TABLE I: Cross sections in fb for e'^e~ annihilation into 

double-charmonium states H1+H2 at -Eboam = 5.3 GeV with- 
out relativistic corrections. The errors are only those from 
variations in the NLO pole mass rric = 1.4 ± 0.2 GeV. Prom 
Ref. [16]. 



H2 \ Hi 


J/V 


i>{2S) 


Vc 


2.31 ± 1.09 


0.96 ± 0.45 




0.96 ± 0.45 


0.40 ± 0.19 


Xco(lP) 


2.28 ± 1.03 


0.95 ± 0.43 



sumed that they are produced by a single virtual pho- 
ton. The four QCD diagrams for the color-singlet pro- 
cess 7* — > cci -I- cci are shown in Fig. 1. We take the 
upper cc pair in Fig. 1 to form a C = — charmonium 
Hi with momentum Pi and the lower cc pair to form 
a, C = + charmonium H2 with momentum P2. There 
are also QED diagrams for 7* cci + cci that can be 
obtained from the QCD diagrams in Fig. 1 by replac- 
ing the virtual gluons by virtual photons, but they are 
suppressed by a factor of a/us- However if one of the 
charmonia is a 1 state, such as a J/'ip, then there are 
the additional QED diagrams in Fig. 2. Although they 
are also suppressed by a factor of a/a^, they are en- 
hanced by a kinematic factor of l/r^, where the variable 
r is defined by 



and therefore can be more important than one might ex- 
pect. 

The a^aj. term in the cross section for e'^e" J/xp+rjc 
was calculated previously by Brodsky and Ji [23]. They 
presented their result in the form of a graph of R ver- 
sus but they did not give an analytic expression 
for the cross section. Our predictions for the double- 
charmonium cross sections without relativistic correc- 
tions are given in Table I. The error bars are those 
associated with the uncertainty in the NLO pole mass 
rric only. The relativistic corrections increase the central 
values of the cross sections by about 2.4 for J/tp + t]c, by 
about 6 for J/ip + rjc{2S) and i^{2S) + rjc, and by about 
13 for ■ip{2S) + ric{2S). Although the total correction fac- 
tor for J/tp + ric is significantly larger than 1, it is the 
product of several modest correction factors that all go 
in the same direction [16]. 

The BELLE Collaboration has recently measured the 
cross section for J/tp + rjc [14]. The J/ip was detected 
through its decays into At"'"/x~ and e+e~, which have a 
combined branching fraction of about 12%. The t]c was 
observed as a peak in the momentum spectrum of the 
J/tjj corresponding to the 2-body process J/tp + r]c- The 
measured cross section is 

a[J/tP + 71,] X B[> 4] = (33t^ ± 9) fb, (3) 

where B [> 4] is the branching fraction for the rjc to decay 
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into at least 4 charged particles. Since B[> 4] < 1, the 
right side of Eq. (3) is a lower bound on the cross section 
for J/t/j + rjc- Updated BELLE cross section is even larger 
than that given in Eq. (3) [22]. 

The lower bound provided by Eq. (3) is about an or- 
der of magnitude larger than the central value 2.3 fb of 
the calculated cross section for J/tp + ric in Table L The 
largest theoretical errors are QCD perturbative correc- 
tions, which we estimate to give an uncertainty of roughly 
60%, the value of rric, which we estimate to give an uncer- 
tainty of roughly 50%, and a relativistic correction that 
we have not been able to quantify with confidence. 

In Ref. [16], we considered various two-charmonium 
processes in which the charmonia have opposite C-parity. 
Complete helicity amplitudes for those processes are also 
presented in Ref. [16]. Liu, He, and Chao also calculated 
the a^a^ terms in the cross sections for e+e^ annihilation 
into J/tp+H, H = r]c, XcO, Xci, and Xc2 [17]. Their results 
are consistent with ours. 

One might suspect that the discrepancy between the- 
ory and experiment exists because of some problem in 
NRQCD factorization approach. Substituting the in- 
put parameters used in above analysis into the formula 
for the rate that is derived from light-front method in 
Ref. [23], we find the the result exactly agrees with that 
from NRQCD [24]. This implies that the discrepancy is 
not from NRQCD, but from pQCD factorization itself. 



III. EXCLUSIVE DOUBLE-J/V PRODUCTION 

Having no conclusive idea for filling the huge gap be- 
tween the data and theory with perturbative corrections, 
we began to question whether the BELLE signal consists 
entirely of J/ip + rjc events. As a possibly missing contri- 
bution, Bodwin, Braaten, and I considered e^e~ annihi- 
lation into double- J/V' states that have the same charge- 
conjugation parity (C-parity) [18, 19]. The strongest mo- 
tivation was the observation that the width of the r]c sig- 
nal measured by the BELLE Collaboration is similar to 
the mass difference between J/il) and r]c- In the BELLE 
fit to the J/ij} momentum distribution, the full width at 
half maximum of the rjc peak is about 0.11 GeV. Since 
the mass difference between the J/^f) and ric is about 0.12 
GeV, there are probably J/^'+J/V' events that contribute 
to the J/i' + rjc signal that is observed by BELLE. 

This process proceeds, at leading order in the QCD 
coupling as, through QED diagrams, shown in Fig. 3, 
that contain two virtual photons. One might expect that 
these cross sections would be much smaller than those 
for charmonia with opposite C-parity because they are 
suppressed by a factor of a^/a^. However, if both char- 
monia have quantum numbers J^'-' = 1 , then there is 
a contribution to the cross section in which each photon 
fragments into a charmonium [18, 19]. The fragmenta- 
tion contribution is enhanced by powers of 1/r [18, 19]. A 
similar, but less dramatic, example is the non-negligible 
QED contribution shown in Fig. 2 for e~^e^ j* ^ 




(c) (d) 



FIG. 3: QED diagrams for the process e+e" —>■ 7*7* cci -I- 
cci. The upper and lower cc pairs evolve into Hi and H2, 
respectively. 

J/ 2jj + rjc process. This enhancement can compensate for 
the suppression factor that is associated with the cou- 
pling constants. 

The photon-fragmentation contributions shown in 
Figs. 3(a) and 3(b) are enhanced because the virtual- 
photon propagators arc of order 1 / instead of order 
l/iJ^gj^jj^. In the amplitude, there arc also two numera- 
tor factors of nic instead of ii'bcam, which arise from the 
cc electromagnetic currents. Hence, the net enhancement 
of the squared amplitude is 1/r^, and the contributions 
to R can be nonzero in the limit r — !■ 0. As r ^ with 
fixed scattering angle 9, the photon-fragmentation con- 
tributions to the cross section factor into the cross section 
for e+e^ "f"f with photon scattering angle 9 and the 
fragmentation probabilities for j ^ Hi and 7 ^ i?2- 
These fragmentation probabilities are nonzero at order 
a only for J^'~' = 1~~ states with hclicitics satisfying 
Ai = — A2 = ±1. The contribution to the ratio R for 
J/tl> + J/tp has the behavior 

ii[J/^(Ai) + J/V(A2)] ~ a^v^, (4) 

where Ai = — A2 = ±1. This asymptotic behavior can be 
compared with that for J/tjj + rjc- 

Rm + Vc] ~ aliv^iry, (5) 

which holds generally for S'-wave final states with oppo- 
site C-parity. In this case, an additional factor arises 
from helicity suppression. The ratio R in Eq. (4) is sup- 
pressed relative to Eq. (5) by a factor of (a/a^)^ « 10~^, 
but the enhancement factor that scales as makes the 
cross sections comparable in magnitude at the energy of 
a B factory. In addition, the differential cross section 
has a sharp peak in both forward and backward regions 
owing to the subproccss e^e~ 77, which results in an- 
other enhancement factor ln(8/r^) in the integrated cross 
section. 

The differential cross section for exclusive double- J/'f/' 
production as a function of a; = cos 9 is shown in Fig. 4. 
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TABLE 11: Cross sections in fb for e^e~ annihilation at 
Sbcam = 5.3 GeV into double-charmonium states H\ + H2 
with C = —1. The errors are only those from variations in the 
pole mass rric — 1.4 ±0.2 GeV. There are additional large er- 
rors associated with perturbative and relativistic corrections, 
as discussed in the text. From Ref. [18, 19]. 



H2 \ Hi 




V'(25) 




8.70 ± 2.94 


7.22 ± 2.44 


V'(25) 




1.50 ± 0.51 




Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0,2 0.4 0.6 0.8 1 

Ixl 



FIG. 4: Differential cross sections da/d\x\ for e^e~ annihila- 
tion at iJbcam = 5.3 GeV into J/tp + J/tp and J/tp + rjc- There 
are large errors associated with perturbative and relativistic 
corrections, as discussed in the text. From Ref. [18, 19]. 

Our predictions for double-charmonium cross sections are 
given in Table II for C = — 1 states. The error bars are 
those associated with the uncertainty in the pole mass nic 
only. The cross sections for the 1~~ states are dominated 
by the photon- fragmentation diagrams in Figs. 3(a) and 
3(b). For rric 1.4 GeV, they contribute 87.5% of the 
J/^+J/ip cross section. The nonfragmentation diagrams 
in Figs. 1(c) and 1(d) contribute 0.7%, while the inter- 
ference term contributes 11.8%. 

The perturbative corrections to the cross section for 
J/tp + J/ip have not yet been calculated. However 
the perturbative corrections to the dominant photon- 
fragmentation diagrams in Figs. 3(a) and 3(b) arc closely 
related to the perturbative correction to the electromag- 
netic annihilation decay rate for J/ip — > e^e^, which 
gives a multiplicative factor 



The perturbative correction to the photon-fragmentation 
terms in the cross section for J /'ip + J /ij) is just the square 
of the expression (6). If we choose the QCD coupling con- 
stant to be Us = 0.25, which corresponds to a renormal- 



ization scale 2mc, then the perturbative correction yields 
a multipHcative factor (0.79)^ = 0.39. The same pertur- 
bative correction factor applies to the cross sections for 
J/ip + ■0(25') and ^{2S) + ^{2S). This perturbative cor- 
rection factor applies only to the leading contributions 
to the cross sections in the limit r ^ 0. However, since 
these contributions are dominant, we conclude that the 
perturbative corrections are likely to decrease the cross 
sections by about a factor of 3. 

The relativistic corrections to the J/ip + J/'0 cross sec- 
tion are significantly smaller than and have the oppo- 
site sign from the relativistic corrections to the J /^p ~\- rjc 
cross section, which are given in Ref. [16]. The relativis- 
tic correction to the fragmentation process is 0.78. For 
TOc = 1.4 GeV, the relativistic corrections to the J/ip + rjc 
cross section are estimated to increase the cross section 
by about a factor 5.5. The large difference in the rela- 
tivistic corrections suggests that there may be large rela- 
tivistic corrections not only to the absolute cross sections 
for double-charmonium production, but also to the ratios 
of those cross sections. 

If we take into account both perturbative and rela- 
tivistic corrections, then the predicted cross section for 
J/ip + J/ip at the B factories is of the same order of 
magnitude as that for J/tJj + rjc- After this proposal was 
made, the BELLE Collaboration looked for the predicted 
double- J/'f/' events. Unfortunately, they did not find such 
events [25]. The non-observation of double- J/i/' events 
strongly suggests there is something wrong in our under- 
standing of factorization in quarkonium production or of 
the relevant production mechanisms. If the factorization 
formalism is valid and we have accounted for the dom- 
inant production mechanisms, then there is no reason 
why J/ip + rjc events should be seen while double- J /ip 
events are not. Note that the two production processes 
depend on essentially the same non-perturbative factor. 
At the moment, we know the relativistic correction seems 
to have an important role in increasing the cross sections 
for J/ijj + ric [16]. The corrections to the two-charmonium 
processes in next-to-leading order in the strong coupling 
constant are still unknown. The NLO result may help us 
with pinning down the origin of the problem. 



IV. EXCLUSIVE J/V-GLUEBALL 
PRODUCTION 

The cross sections for J/ip + rjc, Xco, and ric{2S) re- 
cently measured by the BELLE Collaboration are not 
well understood within pQCD or within NRQCD factor- 
ization. If there is nothing wrong in our factorization 
formula, then the BELLE signal must include something 
else that we have not considered. The non-observation 
of double- J/ V' events strongly suggests there is some- 
thing other than rjc in the signal. This follows from 
the fact that the non-observation of double- J/'0 events 
in the current BELLE data is still consistent with the 
theoretical predictions for the absolute double- J/V' cross 
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FIG. 5: Feynman diagram for 7* ^ i/ + Qj. 



section. On the basis of this observation, Brodsky, Gold- 
haber, and I argued that the signal may include events in 
which there is exclusive J/tp production associated with 
a J^^ = J++ gluebaU Gj with J = 0, 2 [20]. 

Bound states of gluons provide an explicit signature 
of the non-Abelian interactions of quantum chromody- 
namics. In fact, in a model universe without quarks, 
the hadronic spectrum of QCD would consist solely of 
color-singlet glueball states. According to a recent lat- 
tice calculation by Morningstar and Peardon [26], the 
ground-state masses for the J^'~^ — 0++ and 2++ glue- 
balls Gj are 1.73 and 2.40 GeV, respectively. There are 
many excited states in the mass range of charmonium 
spectrum [26]. 

The cross section for e+e^ — > J/ip + Qj might not 
be suppressed compared with cross sections for exclusive 
quarkonium pairs such as 7* — > J -\- rjc that arise from 
the subprocess 7* — > (cc) + (cc). That is because the 
subprocesses 7* — > (cc) -I- (cc) and 7* (cc) + [gg) are 
of the same nominal order in pQCD. Thus, it is possible 
that some portion of the signal observed by BELLE in 
e"'"e~ — > J/V' + X may actually be due to the production 
of J/V) + Gj pairs. 

The amplitude for e+e^ J /tp + Qj at leading twist 
can be expressed as a factorized product of the perturba- 
tive hard-scattering amplitude Tnij* — > QQ+gg) convo- 
luted with the nonperturbative distribution amplitudes 
for the heavy quarkonium and glueball states. A bound 
on the normalization of the distribution amplitude for the 
glueball state can be extracted from a resonance search 
by CUSB in T ^ 7 -I- X [27] following the method used 
in Ref. [28]. 

Our predictions for the upper limits to the exclusive 
charmonium-giueball production cross sections are given 
in Table III. We find that the upper limit to the cross 



section uj/^g^ is comparable to the NRQCD prediction 
of the cross sections for e+e~ J /"ip -\- H for H = rjc 
and Xco, and larger by factor 2 than that for H = r]c{2S), 
suggesting the possibility that a significant fraction of the 
anomalously large cross section measured by BELLE may 
be due to glueballs in association with J/ip production. 

TABLE III: Upper limits to the cross section (Jj/^g^ and 
the ratio (Jj/^g^/aj/^H at y's — 10.6 GeV, assuming that 
Mgg = Mh, where H = rjc, XcO, and rjc{2S). The limits are 
determined by the T — > 7X search of the CUSB Collabora- 
tion [27]. From Ref. [20]. 



Mg„ = Mh 


h = rjc 


XcO 


Vci2S) 


max 


1.4 fb 


1.5 fb 


1.6 fb 


max / 

o"j/^So/o-.//i/''i 


0.63 


0.72 


1.9 



V. DISCUSSION 

The anomalously large cross section for exclusive J/ 4'+ 
rjc production measured by the BELLE Collaboration is 
not well understood within pQCD factorization or within 
the NRQCD factorization formalism. If the BELLE sig- 
nals are made purely of J/ip + rjc events, then either 
pQCD factorization fails or there are important produc- 
tion mechanisms that we have not yet taken into account. 
The former possibility would be a violation of an es- 
tablished pQCD factorization theorem. If the exclusive 
two-charmonium production process is demonstrated to 
violate the factorization theorem, then one must under- 
stand why. Without a fundamental understanding this 
problem, one can not safely use such factorization theo- 
rems, which are a central part of most current particle 
phenomenology 
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